Aplicación de PCR a partir de cultivo en la identificación de subespecies del complejo Mycobacterium Tuberculosis / The use of culture-based PCR for the identification of subspecies of the mycobacterium tuberculosis complex

Los métodos diagnósticos clásicos de tuberculosis (TB) se basan en la utilización de baciloscopía y cultivo. La identificación del agente etiológico desde la positivización del cultivo requiere entre 10 y 15 días, mientras que el empleo de la reacción en cadena de la polimerasa (PCR) disminuye el tiempo a 24 h, lo que permite no solo identificar las subespecies del complejo Mycobacterium tuberculosis (CMTB) sino también diferenciarlas de otras especies ambientales clínicamente importantes (MOTT) facilitando el diagnóstico y tratamiento. El objetivo del presente trabajo fue determinar la utilidad de la PCR en la identificación temprana de las micobacterias pertenecientes al CMTB, a partir de cultivos positivos, de pacientes con sospecha de TB, atendidos en un hospital pediátrico de alta complejidad, durante un período de cuatro años. A cada muestra, se le realizó baciloscopía y cultivo en medio líquido. A los cultivos positivos, una inmunocromatografía lateral (TBIDR) y luego PCR. El 4,6% del total de muestras (510/11.162) pertenecientes a 198 pacientes presentó cultivos positivos. Cuatrocientos veintiseis (84%) correspondieron a muestras respiratorias. El rendimiento de la baciloscopía directa fue del 41% (194/470). Cuatrocientos treinta y ocho (86%) resultaron M. tuberculosis, 21 (4%) Mycobacterium bovis, 7 (1%), M. bovis-BCG y 44 (9%) MOTT. La utilización de medios de cultivos líquidos junto con el empleo de PCR favorecen una rápida orientación microbiológica y constituye una estrategia útil para optimizar el manejo clínico de estas infecciones, desde el punto de vista terapéutico y epidemiológico, especialmente en pediatría (AU)

Classical diagnostic methods for tuberculosis (TB) are based on the use of smear microscopy and culture. The identification of the etiological agent from positive culture requires 10 to 15 days, while the use of the polymerase chain reaction (PCR) reduces the time to 24 h, which allows not only to identify the subspecies of the Mycobacterium tuberculosis complex (MTC) but also to differentiate them from clinically important environmental mycobacteria other than tuberculosis (MOTT), facilitating diagnosis and treatment. The aim of this study was to determine the usefulness of PCR in the early identification of mycobacteria belonging to the MTC, from positive cultures of patients with suspected TB seen in a pediatric tertiary hospital over a 4-year period. For each sample, smear microscopy and culture in liquid medium was performed. Positive cultures were subjected to lateral immunochromatography (TBIDR) and then PCR. Of the total number of samples (510/11,162) belonging to 198 patients, 4.6% showed positive cultures; 426 (84%) were respiratory samples. The direct smear microscopy yield was 41% (194/470). Overall, 438 (86%) were found to be M. tuberculosis, 21 (4%) Mycobacterium bovis, 7 (1%), M. bovis-BCG, and 44 (9%) MOTT. The use of liquid culture media together with the use of PCR favors a rapid microbiological orientation and is a useful strategy to optimize the clinical management of these infections, from a therapeutic and epidemiological point of view, especially in children (AU)

Consistency of Mycobacterium tuberculosis Complex Spoligotyping between the Membrane-Based Method and In Silico Approach.

To tackle the spread of tuberculosis (TB), epidemiological studies are undertaken worldwide to investigate TB transmission chains. Clustered regulatory interspaced short palindromic repeats (CRISPR) locus diversity, also called spoligotyping, is a widely used genotyping assay for the characterization of Mycobacterium tuberculosis complex (MTBC). We compared herein the spoligotyping of MTBC clinical isolates using a membrane-based method (following an initial PCR step) and whole-genome sequencing (WGS)-based method (i.e., in silico spoligotyping). All MTBC strains isolated at the Lyon University Hospital, France, between November 2016 and December 2020 were included (n = 597). Spoligotyping profiles were also used for species identification among the MTBC. Outputs of both methods were analyzed, and discrepant results were investigated thanks to CRISPRbuilder-TB. The overall agreement was 85.7%. Spacer discrepancies observed between the methods were due to the insertion of IS6110 within the direct repeat (DR) sequence upstream or downstream of spacers, mutated DR sequences, or truncated spacers. Discrepancies did not impact species identification. Although spoligotyping-based species identification was inconclusive for 29 isolates, SNP-based phylogeny conducted after WGS allowed the identification of 23 M. tuberculosis (Mtb), 2 M. canettii, and 4 mixed MTBC infections. WGS yielded very few discrepancies compared to membrane-based spoligotyping. Overall agreement was significantly improved (92.4%) by the CRISPR locus reconstruction using CRISPRbuilder-TB for the MTBC isolates with the shared international type 53 in silico spoligotyping. A smooth transition from the membrane-based to the in silico-based genotyping of M. tuberculosis isolates is, therefore, possible for TB diagnosis and epidemiologic survey. IMPORTANCE Whole-genome sequencing (WGS) has profoundly transformed the perspectives of tuberculosis (TB) diagnosis, providing a better discriminatory power to determine relatedness between Mycobacterium tuberculosis complex (MTBC) isolates. Previous genotyping approaches, such as spoligotyping consisting of an initial PCR step followed by reverse dot hybridization, are currently being replaced by WGS. Several pipelines have been developed to extract a spoligotype from WGS data (in silico spoligotyping) allowing for the continuity of MTBC molecular surveys before and after WGS implementation. The present study found very good overall agreement between hybridization to membrane-based spoligotyping and in silico spoligotyping, indicating the possibility of a smooth transition from the traditional to the in silico-based genotyping of MTBC isolates for TB diagnosis and epidemiological survey.

Spoligotyping of Clinical Isolates of Mycobacterium tuberculosis Complex Species in the Oromia Region of Ethiopia.

Background: Tuberculosis (TB) is a leading cause of morbidity and mortality in Ethiopia. Investigation of the Mycobacterium tuberculosis complex (MTBC) species circulating in the Ethiopian population would contribute to the efforts made to control TB in the country. Therefore, this study was conducted to investigate the MTBC species and spoligo patterns in the Oromia region (central) of Ethiopia. Methods: A cross-sectional study design was used to recruit 450 smear positive pulmonary TB (PTB) cases from the Oromia region between September 2017 and August 2018. Mycobacteria were isolated from sputum samples on the Lowenstein Jensen (LJ) medium. Molecular identification of the isolates was performed by spoligotyping. The results of spoligotyping were transferred into a query box in the SITVIT2 database and Run TB-Lineage in the TB Insight website for the identification of spoligo international type (SIT) number and linages of the isolates, respectively. Statistical Product and Service Solutions (SPSS) 20 was applied for statistical analysis. Results: Three hundred and fifteen isolates were grouped under 181 different spoligotype patterns. The most dominantly isolated spoligotype pattern was SIT149 and it consisted of 23 isolates. The majority of the isolates were grouped under Euro-American (EA), East-African-Indian (EAI), and Indo-Oceanic (IO) lineages. These lineages consisted of 79.4, 9.8, and 9.8% of the isolates, respectively. One hundred and sixty-five of the isolates were classified under 31 clustered spoligotypes whereas the remaining 150 were singleton types. Furthermore, 91.1% of the total isolates were classified as orphan types. Clustering of spoligotypes was associated (p < 0.001) with EAI lineage. Conclusion: SIT149 and EA lineage were predominantly isolated from the Oromia region substantiating the findings of the similar studies conducted in other regions of Ethiopia. The observation of significant number of singleton and orphan spoligotypes warrants for additional genetic typing of the isolates using method(s) with a better discriminatory power than spoligotyping.

Impact of pathobiological diversity of Mycobacterium tuberculosis on clinical features and lethal outcome of tuberculosis.

BACKGROUND: Mycobacterium tuberculosis population in Russia is dominated by the notorious Beijing genotype whose major variants are characterized by contrasting resistance and virulence properties. Here we studied how these strain features could impact the progression of pulmonary tuberculosis (TB) concerning clinical manifestation and lethal outcome. RESULTS: The study sample included 548 M. tuberculosis isolates from 548 patients with newly diagnosed pulmonary TB in Omsk, West Siberia, Russia. Strains were subjected to drug susceptibility testing and genotyping to detect lineages, sublineages, and subtypes (within Beijing genotype). The Beijing genotype was detected in 370 (67.5%) of the studied strains. The strongest association with multidrug resistance (MDR) was found for epidemic cluster Beijing B0/W148 (modern sublineage) and two recently discovered MDR clusters 1071-32 and 14717-15 of the ancient Beijing sublineage. The group of patients infected with hypervirulent and highly lethal (in a mouse model) Beijing 14717-15 showed the highest rate of lethal outcome (58.3%) compared to Beijing B0/W148 (31.4%; P = 0.06), Beijing Central Asian/Russian (29.7%, P = 0.037), and non-Beijing (15.2%, P = 0.001). The 14717-15 cluster mostly included isolates from patients with infiltrative but not with fibrous-cavernous and disseminated TB. In contrast, a group infected with low virulent 1071-32-cluster had the highest rate of fibrous-cavernous TB, possibly reflecting the capacity of these strains for prolonged survival and chronicity of the TB process. CONCLUSIONS: The group of patients infected with hypervirulent and highly lethal in murine model 14717-15 cluster had the highest proportion of the lethal outcome (58.3%) compared to the groups infected with Beijing B0/W148 (31.4%) and non-Beijing (15.2%) isolates. This study carried out in the TB high-burden area highlights that not only drug resistance but also strain virulence should be considered in the implementation of personalized TB treatment.

Drug-Resistant Characteristics, Genetic Diversity, and Transmission Dynamics of Rifampicin-Resistant Mycobacterium tuberculosis in Hunan, China, Revealed by Whole-Genome Sequencing.

To gain a deep insight into the additional drug-resistant profiles, genetic diversity, and transmission dynamics of rifampicin-resistant tuberculosis (RR-TB) circulating in Hunan province, drug susceptibility testing and whole-genome-sequencing were performed among RR-TB strains collected from Jan. 2013 to Jun. 2018 in Hunan province. A total of 124 RR-TB strains were recovered successfully and included into the final analysis. Lineage 2.2.1 was the dominant sublineage, accounting for 72.6% (90/124), followed by lineage 4.5 (11.3%, 14/124), lineage 4.4 (8.1%, 10/124), lineage 4.2 (6.5%, 8/124) and lineage 2.2.2 (1.6%, 2/124). Overall, 83.1% (103/124) and 3.2% (4/124) of RR-TB were MDR-TB and XDR-TB, respectively. Nearly 30% of RR-TB isolates were resistant to fluoroquinolones, and 26.6% (33/124) were pre-XDR-TB. Moreover, 30.6% (38/124) of RR-TB strains were identified as phenotypically resistance to pyrazinamide. Totally, 17 clusters containing 48 (38.7%, 48/124) RR-TB strains were identified, ranging in size from 2 to 10 isolates. No significant difference was detected in clustering rate between lineage 2 and lineage 4 (χ2 = 0.027, P = 0.870). Our study revealed the complexity of RR-TB strains circulating in Hunan province with complex additional drug-resistant profile and relatively higher clustering rates. Comprehensive information based on WGS should be used to guide the design of treatment regimens and tailor public interventions. IMPORTANCE Comprehensive information such as genetic background and drug-resistant profile of MTB strains could help to tailor public interventions. However, these data are limited in Hunan province, one of the provinces with high-TB burden in China. So, this study aimed to provide us with deep insight into the molecular epidemiology of RR-TB isolates circulating in Hunan province by combining phenotypic drug susceptibility testing and whole-genome sequencing. To our knowledge, this is the first study to use whole-genome sequencing data of RR-TB strains spanning more than 5 years for molecular epidemiology analysis in Hunan province, which allows us to identify genetic background information and clustered strains more accurately. Our study revealed the complexity of RR-TB strains circulating in Hunan province with complex additional drug-resistant profile and relatively higher clustering rates. Comprehensive information based on WGS should be used to guide the design of treatment regimens and tailor public interventions.

Detection of Mycobacterium tuberculosis multiple strains in sputum samples from patients with pulmonary tuberculosis in south western Uganda using MIRU-VNTR.

Infections with multiple strains of Mycobacterium tuberculosis are now widely recognized as a common occurrence. Identification of patients infected with multiple strains provides both insight into the disease dynamics and the epidemiology of tuberculosis. Analysis of Mycobacterial Interspersed Repetitive Unit-Variable-Number Tandem Repeats (MIRU-VNTR) has been shown to be highly sensitive in detecting multiple M. tuberculosis strains even in sputum. The goal of this study was to identify cases of multiple M. tuberculosis strain infections among patients diagnosed with pulmonary tuberculosis in Southwestern Uganda and assessment of factors associated with multiple strain infections. DNA extracted directly from 78 sputum samples, each from an individual patient, was analyzed using the standard 24 loci MIRU-VNTR typing. Five (6.4%) of the 78 patients were infected with multiple strains of M. tuberculosis with all of them being the newly diagnosed cases while two-thirds of them were co-infected with HIV. Exact regression analysis projected that the natives were more likely to harbor multiple strains (OR; 0.981, 95% CI 0-7.926) as well as those with a high microbial load (OR; 0.390, 95% CI 0-3.8167). Despite these findings being not statistically significant due to the small sample size, this points to a critical component of disease dynamics that has clinical implications and emphasizes a need for a study using a larger cohort. It is also essential to study the potential factors associated with higher risk of exposure to newly diagnosed and HIV positive patients at the community level. In addition, our ability to detect multiple M. tuberculosis strains using the standard 24 loci MIRU-VNTR typing especially with allelic diversity in loci 2059 and 3171, which are excluded from the 15-locus MIRU-VNTR, lead us to recommend the use of this genotyping technique, especially in areas with tuberculosis endemicity similar to this study.

Towards comprehensive understanding of bacterial genetic diversity: large-scale amplifications in Bordetella pertussis and Mycobacterium tuberculosis.

Bacterial genetic diversity is often described solely using base-pair changes despite a wide variety of other mutation types likely being major contributors. Tandem duplication/amplifications are thought to be widespread among bacteria but due to their often-intractable size and instability, comprehensive studies of these mutations are rare. We define a methodology to investigate amplifications in bacterial genomes based on read depth of genome sequence data as a proxy for copy number. We demonstrate the approach with Bordetella pertussis, whose insertion sequence element-rich genome provides extensive scope for amplifications to occur. Analysis of data for 2430 B. pertussis isolates identified 272 putative amplifications, of which 94 % were located at 11 hotspot loci. We demonstrate limited phylogenetic connection for the occurrence of amplifications, suggesting unstable and sporadic characteristics. Genome instability was further described in vitro using long-read sequencing via the Nanopore platform, which revealed that clonally derived laboratory cultures produced heterogenous populations rapidly. We extended this research to analyse a population of 1000 isolates of another important pathogen, Mycobacterium tuberculosis. We found 590 amplifications in M. tuberculosis, and like B. pertussis, these occurred primarily at hotspots. Genes amplified in B. pertussis include those involved in motility and respiration, whilst in M. tuberuclosis, functions included intracellular growth and regulation of virulence. Using publicly available short-read data we predicted previously unrecognized, large amplifications in B. pertussis and M. tuberculosis. This reveals the unrecognized and dynamic genetic diversity of B. pertussis and M. tuberculosis, highlighting the need for a more holistic understanding of bacterial genetics.

Application of Amplicon-Based Targeted NGS Technology for Diagnosis of Drug-Resistant Tuberculosis Using FFPE Specimens.

Next-generation sequencing (NGS) enables rapid identification of common and rare drug-resistant genetic variations from tuberculosis (TB) patients' sputum samples and MTB isolates. However, whether this technology is effective for formalin-fixed and paraffin-embedded (FFPE) tissues remains unclear. An amplicon-based targeted NGS sequencing panel was developed to predict susceptibility to 9 antituberculosis drugs, including 3 first-line drugs, by directly detecting FFPE tissues. A total of 178 tissue samples from TB patients who underwent phenotypic drug susceptibility test were retrospectively tested from January 2017 to October 2019 in the Department of Pathology, Beijing Chest Hospital, China. Phenotypic drug susceptibility test results were used as the reference standard. We identified 22 high-quality mutations from 178 FFPE tissue samples, including 15 high+moderate+minimal confidence-level mutations associated with drug resistance (rpoB D435V, S450F/L; KatG S315T; inhA-fabG promoter c-15t; embB G406S, M306V; rpsL K43R, K88R, rrs a1401g, a514c; gyrA D94G/Y/A, A90V), 6 mutations not associated with resistance (rpoB D435Y, H445S, L430P, L452P; embB G406A/D), and one mutation site embB M306I defined as indeterminate. Compared to the phenotypic method, sensitivities (95% CI) for rifampicin, isoniazid, and ethambutol were 96% (79.65-99.90%), 93.55% (78.58-99.21%), and 71.43% (35.24-92.44%), respectively; while for second-line drugs, it varied from 23.53% (9.05-47.77%) for capreomycin to 86.84% (72.20-94.72%) for streptomycin. Specificities for all drugs were satisfactory (>94.51%). Therefore, important pathological FFPE tissue samples, despite partially degraded DNA, can be used as essential specimens for molecular diagnosis of drug resistant TB by amplicon-based targeted NGS technology. IMPORTANCE Amplicon-based targeted NGS technology focuses on a set of gene mutations of known or suspected associations with drug susceptibility in Mycobacterium tuberculosis (MTB). This method offers many benefits, such as low sequencing cost, easy customization, high throughput, shorter testing time and not culture dependent. Formalin-fixed and paraffin-embedded (FFPE) tissues are important pathological specimen in diagnosing tuberculous disease because they are noninfectious and provide excellent preservation of tissue morphology with low storage cost. However, the performance of amplicon-based targeted NGS method on FFPE samples has not been reported yet. Therefore, we evaluated the performance of this method using FFPE samples collected from January 2017 to October 2019 in the Department of Pathology, Beijing Chest Hospital, China. We demonstrate that the amplicon-based targeted NGS method performs excellent on FFPE samples, and it can be applied to pathological diagnosis of drug resistant tuberculosis.

Whole-genome sequencing of Mycobacterium tuberculosis for prediction of drug resistance.

Whole-genome sequencing (WGS) has shown tremendous potential in rapid diagnosis of drug-resistant tuberculosis (TB). In the current study, we performed WGS on drug-resistant Mycobacterium tuberculosis isolates obtained from Shanghai (n = 137) and Russia (n = 78). We aimed to characterise the underlying and high-frequency novel drug-resistance-conferring mutations, and also create valuable combinations of resistance mutations with high predictive sensitivity to predict multidrug- and extensively drug-resistant tuberculosis (MDR/XDR-TB) phenotype using a bootstrap method. Most strains belonged to L2.2, L4.2, L4.4, L4.5 and L4.8 lineages. We found that WGS could predict 82.07% of phenotypically drug-resistant domestic strains. The prediction sensitivity for rifampicin (RIF), isoniazid (INH), ethambutol (EMB), streptomycin (STR), ofloxacin (OFL), amikacin (AMK) and capreomycin (CAP) was 79.71%, 86.30%, 76.47%, 88.37%, 83.33%, 70.00% and 70.00%, respectively. The mutation combination with the highest sensitivity for MDR prediction was rpoB S450L + rpoB H445A/P + katG S315T + inhA I21T + inhA S94A, with a sensitivity of 92.17% (0.8615, 0.9646), and the mutation combination with highest sensitivity for XDR prediction was rpoB S450L + katG S315T + gyrA D94G + rrs A1401G, with a sensitivity of 92.86% (0.8158, 0.9796). The molecular information presented here will be of particular value for the rapid clinical detection of MDR- and XDR-TB isolates through laboratory diagnosis.

Accuracy of an amplicon-sequencing nanopore approach to identify variants in tuberculosis drug-resistance-associated genes.

A rapid and accurate diagnostic assay represents an important means to detect Mycobacterium tuberculosis, identify drug-resistant strains and ensure treatment success. Currently employed techniques to diagnose drug-resistant tuberculosis include slow phenotypic tests or more rapid molecular assays that evaluate a limited range of drugs. Whole-genome-sequencing-based approaches can detect known drug-resistance-conferring mutations and novel variations; however, the dependence on growing samples in culture, and the associated delays in achieving results, represents a significant limitation. As an alternative, targeted sequencing strategies can be directly performed on clinical samples at high throughput. This study proposes a targeted sequencing assay to rapidly detect drug-resistant strains of M. tuberculosis using the Nanopore MinION sequencing platform. We designed a single-tube assay that targets nine genes associated with drug resistance to seven drugs and two phylogenetic-determining regions to determine strain lineage and tested it in nine clinical isolates and six sputa. The study's main aim is to calibrate MinNION variant calling to detect drug-resistance-associated mutations with different frequencies to match the accuracy of Illumina (the current gold-standard sequencing technology) from both culture and sputum samples. After calibrating Nanopore MinION variant calling, we demonstrated 100% agreement between Illumina WGS and our MinION set up to detect known drug resistance and phylogenetic variants in our dataset. Importantly, other variants in the amplicons are also detected, decreasing the recall. We identify minority variants and insertions/deletions as crucial bioinformatics challenges to fully reproduce Illumina WGS results.

Strong Increase in Moxifloxacin Resistance Rate among Multidrug-Resistant Mycobacterium tuberculosis Isolates in China, 2007 to 2013.

We designed this study to determine the trend of moxifloxacin resistance among multidrug-resistant tuberculosis (MDR-TB) patients from 2007 to 2013 in China to inform the composition of multidrug-resistant/rifampicin-resistant tuberculosis (MDR/RR-TB) treatment regimens. We assessed moxifloxacin resistance among MDR-TB isolates collected in national drug resistance surveys in 2007 and 2013 that included 3,634 smear-positive and 7,206 culture-positive pulmonary tuberculosis patients, respectively. Moxifloxacin susceptibility was examined by a Mycobacterium growth indicator tube (MGIT) 960 for the 2007 isolates, and by the minimum inhibitory concentration (MIC) method for the 2013 isolates, at both breakpoints 0.5 and 2.0 µg/mL. Risk factors were explored through multivariable log-binominal regression analysis. Mutations in gyrA and gyrB for part of the isolates were also studied through sequencing. Of 401 MDR strains isolated in 2007, moxifiloxacin resistance could be determined for 319 (79.6%): 41 (12.9%) and 10 (3.1%) were resistant at 0.5 and 2.0 µg/mL, respectively. Of 365 MDR strains isolated in 2013, 338 (92.6%) could be analyzed: 140 (41.4%) and 79 (23.4%) were resistant at 0.5 and 2.0 µg/mL. For patients in 2007, no characteristics were significantly associated with moxifloxacin resistance. For patients in 2013, patients aged ≥60 years (adjusted prevalence ratio [aPR], 1.46; 95% confidence interval [CI], 1.10 to 1.93) were more likely to have resistance at 0.5 µg/mL, whereas those residing in eastern China compared to those in central China had an increased risk of resistance at both 0.5 (aPR, 1.85; 95% CI, 1.38 to 2.48) and 2.0 µg/mL (aPR, 2.14; 95% CI, 1.35 to 3.40). Sequencing results were obtained for 245 and 266 MDR-TB isolates in 2007 and 2013, respectively. In total, 34 of 38 (89.5%) and 89 of 104 (85.6%) of 2007 and 2013 moxifloxacin-resistant (0.5 µg/mL) MDR-TB strains had mutations in the gyrA and gyrB gene, respectively. Asp94Gly was the most common mutation among 2007 (11 of 38, 28.9%) and 2013 isolates (24 of 104, 23.1%) and conferred high-level moxifloxacin resistance. Moxifloxacin resistance among MDR-TB patients in China increased from modest to high from 2007 to 2013. Moxifloxacin should be used carefully as a potentially effective drug for composing MDR/RR-TB regimens especially for elderly patients in China. Individual susceptibility testing especially rapid molecular-based assays should be conducted to confirm the susceptibility to moxifloxacin. IMPORTANCE China is one of the high-burden countries for multidrug-resistant/rifampicin-resistant tuberculosis (MDR/RR-TB). Moxifloxacin is one of the critical antituberculosis drugs for MDR/RR-TB treatment. Susceptibility to moxifloxacin is therefore very important to compose effective regimens and to provide protection against development of resistance of companion drugs such as bedaquiline and linezolid. There are, however, no nationally representative data on moxifloxacin resistance among MDR/RR-TB cases in China. Therefore, we assessed the resistance prevalence for moxifloxacin among MDR-TB strains isolated in national drug resistance surveys in 2007 and 2013 that covered 72 sites around the country. We demonstrate that the prevalence of moxifloxacin resistance in MDR-TB isolates increased from modest to high, which should prompt the national tuberculosis program to use moxifloxacin cautiously in second-line regimens to treat MDR/RR-TB unless susceptibility can be laboratory-confirmed.

Transcriptomic Characterization of Tuberculous Sputum Reveals a Host Warburg Effect and Microbial Cholesterol Catabolism.

The crucial transmission phase of tuberculosis (TB) relies on infectious sputum and yet cannot easily be modeled. We applied one-step RNA sequencing (RNA-Seq) to sputum from infectious TB patients to investigate the host and microbial environments underlying transmission of Mycobacterium tuberculosis. In such TB sputa, compared to non-TB controls, transcriptional upregulation of inflammatory responses, including an interferon-driven proinflammatory response and a metabolic shift toward glycolysis, was observed in the host. Among all bacterial sequences in the sputum, approximately 1.5% originated from M. tuberculosis, and its transcript abundance was lower in HIV-1-coinfected patients. Commensal bacterial abundance was reduced in the presence of M. tuberculosis infection. Direct alignment to the genomes of the predominant microbiota species also reveals differential adaptation, whereby firmicutes (e.g., streptococci) displayed a nonreplicating phenotype with reduced transcription of ribosomal proteins and reduced activities of ATP synthases, while Neisseria and Prevotella spp. were less affected. The transcriptome of sputum M. tuberculosis more closely resembled aerobic replication and shared similarity in carbon metabolism to in vitro and in vivo models with significant upregulation of genes associated with cholesterol metabolism and downstream propionate detoxification pathways. In addition, and counter to previous reports on intracellular M. tuberculosis infection in vitro, M. tuberculosis in sputum was zinc, but not iron, deprived, and the phoP loci were also significantly downregulated, suggesting that the pathogen is likely extracellular in location. IMPORTANCE Although a few studies have described the microbiome composition of TB sputa based on 16S ribosomal DNA, these studies did not compare to non-TB samples and the nature of the method does not allow any functional inference. This is the first study to apply such technology using clinical specimens and obtained functional transcriptional data on all three aspects simultaneously. We anticipate that an improved understanding on the biological interactions in the respiratory tract may also allow novel interventions, such as those involving microbiome manipulation or inhibitor targeting disease-specific metabolic pathways.

SNPPar: identifying convergent evolution and other homoplasies from microbial whole-genome alignments.

Homoplasic SNPs are considered important signatures of strong (positive) selective pressure, and hence of adaptive evolution for clinically relevant traits such as antibiotic resistance and virulence. Here we present a new tool, SNPPar, for efficient detection and analysis of homoplasic SNPs from large whole genome sequencing datasets (>1000 isolates and/or >100 000 SNPs). SNPPar takes as input an SNP alignment, tree and annotated reference genome, and uses a combination of simple monophyly tests and ancestral state reconstruction (ASR, via TreeTime) to assign mutation events to branches and identify homoplasies. Mutations are annotated at the level of codon and gene, to facilitate analysis of convergent evolution. Testing on simulated data (120 Mycobacterium tuberculosis alignments representing local and global samples) showed SNPPar can detect homoplasic SNPs with very high specificity (zero false-positives in all tests) and high sensitivity (zero false-negatives in 89 % of tests). SNPPar analysis of three empirically sampled datasets (Elizabethkingia anophelis, Burkholderia dolosa and M. tuberculosis) produced results that were in concordance with previous studies, in terms of both individual homoplasies and evidence of convergence at the codon and gene levels. SNPPar analysis of a simulated alignment of ~64 000 genome-wide SNPs from 2000 M. tuberculosis genomes took ~23 min and ~2.6 GB of RAM to generate complete annotated results on a laptop. This analysis required ASR be conducted for only 1.25 % of SNPs, and the ASR step took ~23 s and 0.4 GB of RAM. SNPPar automates the detection and annotation of homoplasic SNPs efficiently and accurately from large SNP alignments. As demonstrated by the examples included here, this information can be readily used to explore the role of homoplasy in parallel and/or convergent evolution at the level of nucleotide, codon and/or gene.

A Rapid Drug Resistance Genotyping Workflow for Mycobacterium tuberculosis, Using Targeted Isothermal Amplification and Nanopore Sequencing.

Phenotypic drug susceptibility testing (DST) for tuberculosis (TB) requires weeks to yield results. Although molecular tests rapidly detect drug resistance-associated mutations (DRMs), they are not scalable to cover the full genome and the many DRMs that can predict resistance. Whole-genome sequencing (WGS) methods are scalable, but if conducted directly on sputum, typically require a target enrichment step, such as nucleic acid amplification. We developed a targeted isothermal amplification-nanopore sequencing workflow for rapid prediction of drug resistance of TB isolates. We used recombinase polymerase amplification (RPA) to perform targeted isothermal amplification (37°C for 90 min) of three regions within the Mycobacterium tuberculosis genome, followed by nanopore sequencing on the MinION. We tested 29 mycobacterial genomic DNA extracts from patients with drug-resistant (DR) TB and compared our results to those of WGS by Illumina and phenotypic DST to evaluate the accuracy of prediction of resistance to rifampin and isoniazid. Amplification by RPA showed fidelity equivalent to that of high-fidelity PCR (100% concordance). Nanopore sequencing generated DRM predictions identical to those of WGS, with considerably faster sequencing run times of minutes rather than days. The sensitivity and specificity of rifampin resistance prediction for our workflow were 96.3% (95% confidence interval [CI], 81.0 to 99.9%) and 100.0% (95% CI, 15.8 to 100.0%), respectively. For isoniazid resistance prediction, the sensitivity and specificity were 100.0% (95% CI, 86.3 to 100.0%) and 100.0% (95% CI, 39.8 to 100.0%), respectively. The workflow consumable costs per sample are less than £100. Our rapid and low-cost drug resistance genotyping workflow provides accurate prediction of rifampin and isoniazid resistance, making it appropriate for use in resource-limited settings. IMPORTANCE Current methods for diagnosing drug-resistant tuberculosis are time consuming, resulting in delays in patients receiving treatment and in transmission onwards. They also require a high level of laboratory infrastructure, which is often only available at centralized facilities, resulting in further delays to diagnosis and additional barriers to deployment in resource-limited settings. This article describes a new workflow that can diagnose drug-resistant TB in a shorter time, with less equipment, and for a lower price than current methods. The amount of TB DNA is first increased without the need for bulky and costly thermocycling equipment. The DNA is then read using a portable sequencer called a MinION, which indicates whether there are tell-tale changes in the DNA that indicate whether the TB strain is drug resistant. Our workflow could play an important role in the future in the fight against the public health challenge that is TB drug resistance.

Comparative transcriptomic analysis of THP-1-derived macrophages infected with Mycobacterium tuberculosis H37Rv, H37Ra and BCG.

Tuberculosis (TB) remains a worldwide healthcare concern, and the exploration of the host-pathogen interaction is essential to develop therapeutic modalities and strategies to control Mycobacterium tuberculosis (M.tb). In this study, RNA sequencing (transcriptome sequencing) was employed to investigate the global transcriptome changes in the macrophages during the different strains of M.tb infection. THP-1 cells derived from macrophages were exposed to the virulent M.tb strain H37Rv (Rv) or the avirulent M.tb strain H37Ra (Ra), and the M.tb BCG vaccine strain was used as a control. The cDNA libraries were prepared from M.tb-infected macrophages and then sequenced. To assess the transcriptional differences between the expressed genes, the bioinformatics analysis was performed using a standard pipeline of quality control, reference mapping, differential expression analysis, protein-protein interaction (PPI) networks, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Q-PCR and Western blot assays were also performed to validate the data. Our findings indicated that, when compared to BCG or M.tb H37Ra infection, the transcriptome analysis identified 66 differentially expressed genes in the M.tb H37Rv-infected macrophages, out of which 36 genes were up-regulated, and 30 genes were down-regulated. The up-regulated genes were associated with immune response regulation, chemokine secretion, and leucocyte chemotaxis. In contrast, the down-regulated genes were associated with amino acid biosynthetic and energy metabolism, connective tissue development and extracellular matrix organization. The Q-PCR and Western blot assays confirmed increased expression of pro-inflammatory factors, altered energy metabolic processes, enhanced activation of pro-inflammatory signalling pathways and increased pyroptosis in H37Rv-infected macrophage. Overall, our RNA sequencing-based transcriptome study successfully identified a comprehensive, in-depth gene expression/regulation profile in M.tb-infected macrophages. The results demonstrated that virulent M.tb strain H37Rv infection triggers a more severe inflammatory immune response associated with increased tissue damage, which helps in understanding the host-pathogen interaction dynamics and pathogenesis features in different strains of M.tb infection.

Mycobacterium tuberculosis strain lineage in mixed tribal population across India and Andaman Nicobar Island.

In India, the tribal population constitutes almost 8.6% of the nation's total population. Despite their large presence, there are only a few reports available on Mycobacterium tuberculosis (M. tb) strain prevalence in Indian tribal communities considering the mobile nature of this population and also the influence of the mainstream populations they coexist within many areas for their livelihood. This study attempts to provide critical information pertaining to the TB strain diversity, its public health implications, and distribution among the tribal population in eleven Indian states and Andaman & Nicobar (A&N) Island. The study employed a population-based molecular approach. Clinical isolates were received from 66 villages (10 states and Island) and these villages were selected by implying situation analysis. A total of 78 M. tb clinical isolates were received from 10 different states and A&N Island. Among these, 16 different strains were observed by spoligotyping technique. The major M. tb strains spoligotype belong to the Beijing, CAS1_DELHI, and EAI5 family of M. tb strains followed by EAI1_SOM, EAI6_BGD1, LAM3, LAM6, LAM9, T1, T2, U strains. Drug-susceptibility testing (DST) results showed almost 15.4% of clinical isolates found to be resistant to isoniazid (INH) or rifampicin (RMP) + INH. Predominant multidrug-resistant (MDR-TB) isolates seem to be Beijing strain. Beijing, CAS1_DELHI, EAI3_IND, and EAI5 were the principal strains infecting mixed tribal populations across India. Despite the small sample size, this study has demonstrated higher diversity among the TB strains with significant MDR-TB findings. Prevalence of Beijing MDR-TB strains in Central, Southern, Eastern India and A&N Island indicates the transmission of the TB strains.

A bioinformatics pipeline for Mycobacterium tuberculosis sequencing that cleans contaminant reads from sputum samples.

Next-Generation Sequencing (NGS) is widely used to investigate genomic variation. In several studies, the genetic variation of Mycobacterium tuberculosis has been analyzed in sputum samples without previous culture, using target enrichment methodologies for NGS. Alignments obtained by different programs generally map the sequences under default parameters, and from these results, it is assumed that only Mycobacterium reads will be obtained. However, variants of interest microorganism in clinical samples can be confused with a vast collection of reads from other bacteria, viruses, and human DNA. Currently, there are no standardized pipelines, and the cleaning success is never verified since there is a lack of rigorous controls to identify and remove reads from other sputum-microorganisms genetically similar to M. tuberculosis. Therefore, we designed a bioinformatic pipeline to process NGS data from sputum samples, including several filters and quality control points to identify and eliminate non-M. tuberculosis reads to obtain a reliable genetic variant report. Our proposal uses the SURPI software as a taxonomic classifier to filter input sequences and perform a mapping that provides the highest percentage of Mycobacterium reads, minimizing the reads from other microorganisms. We then use the filtered sequences to perform variant calling with the GATK software, ensuring the mapping quality, realignment, recalibration, hard-filtering, and post-filter to increase the reliability of the reported variants. Using default mapping parameters, we identified reads of contaminant bacteria, such as Streptococcus, Rhotia, Actinomyces, and Veillonella. Our final mapping strategy allowed a sequence identity of 97.8% between the input reads and the whole M. tuberculosis reference genome H37Rv using a genomic edit distance of three, thus removing 98.8% of the off-target sequences with a Mycobacterium reads loss of 1.7%. Finally, more than 200 unreliable genetic variants were removed during the variant calling, increasing the report's reliability.

Population structure, biogeography and transmissibility of Mycobacterium tuberculosis.

Mycobacterium tuberculosis is a clonal pathogen proposed to have co-evolved with its human host for millennia, yet our understanding of its genomic diversity and biogeography remains incomplete. Here we use a combination of phylogenetics and dimensionality reduction to reevaluate the population structure of M. tuberculosis, providing an in-depth analysis of the ancient Indo-Oceanic Lineage 1 and the modern Central Asian Lineage 3, and expanding our understanding of Lineages 2 and 4. We assess sub-lineages using genomic sequences from 4939 pan-susceptible strains, and find 30 new genetically distinct clades that we validate in a dataset of 4645 independent isolates. We find a consistent geographically restricted or unrestricted pattern for 20 groups, including three groups of Lineage 1. The distribution of terminal branch lengths across the M. tuberculosis phylogeny supports the hypothesis of a higher transmissibility of Lineages 2 and 4, in comparison with Lineages 3 and 1, on a global scale. We define an expanded barcode of 95 single nucleotide substitutions that allows rapid identification of 69 M. tuberculosis sub-lineages and 26 additional internal groups. Our results paint a higher resolution picture of the M. tuberculosis phylogeny and biogeography.

Drug resistant tuberculosis cases from the Copperbelt province and Northern regions of Zambia: Genetic diversity, demographic and clinical characteristics.

Tuberculosis (TB) caused by Mycobacterium tuberculosis remains a major cause of death worldwide. Diverse genotypes have been demonstrated to drive the epidemiology of drug resistant (DR-) TB globally. Currently, there is limited knowledge on the genotypes and transmission dynamics of M. tuberculosis in Zambia. This study aimed to describe the genotypes of DR-TB from the Copperbelt and Northern regions of Zambia. Molecular typing tools of insertion sequence 6110-restriction fragment length polymorphism (IS6110-RFLP) and spacer oligonucleotide typing (spoligotyping) were applied. We demonstrate that diverse genotypes are associated with DR-TB in Zambia. The predominant genotype was lineage 4; other strains belonged to lineage 2 and 3. Genotypes previously identified as driving the epidemiology of drug susceptible TB have been identified as drivers of DR-TB. Genotyping analysis showed clustering of strains among patients from different regions of the country; suggesting that DR-TB is widespread. Molecular findings combined with phenotypic and epidemiologic findings play a critical role in identifying circulating genotypes and possible transmission chains. Clustering of drug resistant strains was demonstrated to be 48% and 86% according to IS6110-RFLP and spoligotyping, respectively. However, gaps in clinical and demographic data skew the interpretation, and call for data collection policy improvements.

Whole genome sequencing of clinical samples reveals extensively drug resistant tuberculosis (XDR TB) strains from the Beijing lineage in Nigeria, West Africa.

Multi-drug (MDR) and extensively drug-resistant (XDR) tuberculosis (TB) continues to be a global public health problem especially in high TB burden countries like Nigeria. Many of these cases are undetected and go on to infect high risk individuals. Clinical samples from positive rifampicin resistant Xpert®MTB/Rif assay were subjected to direct whole genome sequencing and bioinformatics analysis to identify the full antibiotics resistance and lineage profile. We report two (2) XDR TB samples also belonging to the East-Asian/Beijing family of lineage 2 Mycobacterium tuberculosis complex from clinical samples in Nigeria. Our findings further reveal the presence of mutations that confer resistance to first-line drugs (rifampicin, isoniazid, ethambutol and pyrazanimide), second-line injectables (capreomycin, streptomycin, kanamycin and/or amikacin) and at least one of the fluoroquinolones (ofloxacin, moxifloxacin, levofloxacin and/or ciprofloxacin) in both samples. The genomic sequence data from this study not only provide the first evidence of XDR TB in Nigeria and West Africa, but also emphasize the importance of WGS in accurately detecting MDR and XDR TB, to ensure adequate and proper management treatment regimens for affected individuals. This will greatly aid in preventing the spread of drug resistance TB in high burden countries.